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I. INTRODUCTION 



preserving the gap in the excitation spectrum. 



One of the more intriguing properties of the cuprate su- 
perconductors is the existence of a "pseudogap" regime in 
the normal phase, which develops when these materials 
are underdoped and at temperatures below a characteris- 
tic temperature T* (see Fig. 1). Numerous experiments, 
including NMR,El .neutron scattering,! infrared and op- 
tical conductivity!!! transport,!] Raman scatteringjlspe- 
cific heat,El angle-resolved photoemission. (ARPES), El and 
scanning tunneling spectroscopy (STS)|j indicate that in 
this regime there is a gap in the excitation spectrum, 
while transport measurements show that the sample is 
normal and that the superfluid density is zero.El The 
ARPES and STS experiments on B^S^CaCu^Os+a have 
established that the pseudogap evolves smoothly from 
the superconducting gap. More significantly, the ARPES 
results show that the pseudogap in the normal phase pos- 
sesses the same d x 2_ y 2 symmetry as the superconducting 
gap; the node along the (tt/2, tt/2) direction in the super- 
conducting phase evolves into an extended gapless region 
with increasing temperature in the pseudogap regime, 
while the gap maximum in the (-7T, 0) direction has only a 
very weak temperature dependence. Taken collectively, 
these results demonstrate that the pseudogap in the nor- 
mal phase is a remnant of the quasiparticle gap in the 
superconducting phase. 

The theoretical challenge in understanding the pseudo- 
gap regime is then to reconcile the existence of a quasi- 
particle gap with the absence of superconducting order. 
The usual strategy in studying superconductivity is to 
start from a normal phase consisting of well-defined Lan- 
dau quasiparticles and to then ask what interactions lead 
to a superconducting instability. Here, the experiments 
suggest that we start from the superconducting phase — 
skirting the issue of the origin of the pairing — and ask 
what mechanisms would destroy phase coherence while 
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FIG. 1. Schematic phase diagram of a cuprate supercon- 
ductor as a function of temperature T and doping x. AFM 
denotes the antiferromagnetic insulating phase. 

The mechanism which we explore in this,paper, orig- 
inally suggested by Emery and Kivelson|L3 is that the 
local superconducting order parameter A(x) = Aoe l9 ^ 
has a well defined amplitude A (set by energies of or- 
der T*) but a fluctuating phase 9, such that (A(x)) = 
A (e* e( - x - ) ) = in the pseudogap regime. The tempera- 
ture scale for phase fluctuations T phasc is controlled by 
the zero-temperature superfluid stiffness n s (0)/4m; the 
actual critical temperature T c = min{T phasc , T MF }, with 
T MF the mean-field transition temperature predicted in 
BCS theory. Emery and Kivclson argued that in con- 
ventional low-T c superconductors, which have a large su- 
perfluid stiffness, T phaso > T c , so that phase fluctua- 
tions are irrelevant and the transition is BCS-like. How- 
ever, cuprate superconductors are doped Mott insulators 
which have a small superfluid stiffness (of order the dop- 
ing x), so that T phasc ~ T c and the transition is con- 
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trolled by phase-ordering. This naturally explains the 
Uemura scalingjlil T c ~ n s (0)/m, observed in the cuprate 
superconductors. In addition, the reduction of the phase- 
stiffness n s (T)/Am is approximately linear in tempera- 
ture due to the d-wave symmetry, almost all the way 
to T c ,t3 in contrast to conventional s-wave supercondrie* 
tors, where n s (T)/4m rapidly drops to zero near T c .t3 
Therefore the temperature range for which n s (T) / Am is 
a small number is much wider than in the conventional 
superconductors, and the effect of phase fluctuations can 
be observed even far from T c . 

In this paper we investigate the effects of phase fluc- 
tuations on the fermion single-particle properties and 
the thermodynamics near the superconducting transi- 
tion, without regard to the details of a microscopic model. 
In addition to the smallness of n s (T) / Am, we assume 
that the superconducting gap and the normal state pseu- 
dogap have a d-wave symmetry and that the system is 
two-dimensional due to the weakness of the interlayer 
coupling. We find that phase fluctuations, and in partic- 
ular transverse phase fluctuations, can have a pronounced 
effect on the spectral function, density of states, and spe- 
cific heat near the transition. 

A brief overview of the paper is as follows: In Sec. II we 
construct a low-energy effective theory of fermion quasi- 
particles coupled to the classical phase fluctuations by 
integrating out the fast momentum degrees of freedom. 
In Sec. Ill we obtain the self-energy of the (neutral) 
fermion quasiparticles and calculate the single-particle 
density of states and the spectral broadening. We ob- 
serve that the peak in the spectral function and the den- 
sity of states near the gap is significantly broadened by 
the phase fluctuations (vortex pair unbinding) near the 
transition temperature, exhibiting a pseudogap-like be- 
havior. In Sec. IV we consider the effect of the phase fluc- 
tuations on the specific heat near the transition. We find 
that the vortex-antivortex unbinding transition causes a 
visible peak in the specific heat just above the transi- 
tion temperature due to the short coherence length of the 
materials. Appendix A is a compilation of some useful 
results on the Berezinskii-Kosterlitz-Thouless transition, 
and Appendix B contains some technical details on the 
calculation of the fermion self-energy. 
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with the interaction strength g > 0. The pairing field 
A(x, t) = |A(x, r)|e ie ( x ' T '; in what follows we will as- 
sume that the amplitude of the order parameter is con- 
stant, |A(x, t) = A, and concentrate on the phase de- 
gree of freedom, 8. Also we have not considered the effect 
of Coulomb interactions. In order to couple the fcrmions 
to the 9 field more explicitly, we perform a gauge transfor- 
mation Vv(x, t) = Co.(x, T)e~ i p?' T ^ 2 , with tp a the field 
operators for neutral fermionsEj Then we obtain a new 
action for the neutral fermions as 
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The last term in Eq. (2.3) is a constant and will be 
dropped in what follows. 




II. EFFECTIVE LOW-ENERGY THEORY 

In this section we derive a low-energy effective theory 
of fermion quasiparticles coupled to the phase fluctua- 
tions of the superconducting order parameter, assuming 
that the amplitude fluctuations are negligible near T c . 
For simplicity, we consider an s-wave symmetry of the 
gap in this section. Extensions to other gap symmetries 
are straightforward. 

We begin with the mean-field BCS model whose par- 
tition function is given by 



FIG. 2. Hilbert space of the low-energy effective theory. 
After the inner and outer side of the thin momentum shell 
degree of freedom has been integrated out, we are left with 
a theory of fermion quasiparticles living in a thin momentum 
shell of thickness 2A. The dashed line denotes the Fermi sur- 
face. The direction of the Fermi velocity v^(k) is in parallel 
with the momentum k of the fermion quasiparticle ip(k). 



As it stands, Eq. (2.3) is hard to study, so we integrate 
out the fast momentum degrees of freedom and focus on 
the low-energy properties of the system.Ej Namely, ip — > 
ip + ip' where ^'(x) = X)||k|-fe F |>A V^e* 1 *'* is the fast mo- 
mentum degree of freedom which is to be integrated out 
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of the partition function and i/)(x) = S||k|-fc F |<A x 
is the low-energy quasiparticle. Then 

= Af J VipV^V8e- s ° !t 7 (2.4) 

and we are left with the momentum degrees of freedom in 
the vicinity of the Fermi surface if we take A fci?. Since 
we will be considering the effect of fluctuating vortices on 
the fermion self-energy, and the characteristic vortex core 
size is of order the coherence length £o, we will take A of 
order l/£o- We are therefore working in the limiLij F £ » 
1; in the cuprates it is estimated that kp£,o ~ 10ji3 so our 
approximation is appropriate for these materials. 

The effective action of the fermions near the Fermi 
surface and the phase field is 

&ff = Y1 [r5^Vi(k,w)^(k,w)(iw-fc a /2m + /i) 

tr k.cj 

+ J d 2 x dTipt(x, T)(id T 6/2 
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+ T^A^{(k,w)V>{(-k, -w) + h.c. (2.5) 



Here we have neglected higher orders in q/k <C 1 since 
k w kF and |q| < A, and we have taken cj.c<j w n/, 
the total number of fermions. Since k is close to the 
Fermi surface, we may linearize the fermionic spectrum 
as k 2 /2m — \i = £k ~ w F(|k — kp). The resulting ef- 
fective action is expressed in terms of the Nambu spinor 
notation, ip = (ip^ , ip j ) , as 
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with Go the bare Green's function for the neutral 
fermions, 
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The microscopic physics relevant to the length scale 
L < £o (the vortex core structure, for instance) has been 
integrated out to renormalize such quantities as the ef- 
fective mass m and the residual interactions between the 
quasiparticles which we have not included here. There- 
fore, in this effective theory the vortices are treated as 
point-like entities. 

In studying the finite temperature superconductor to 
normal metal transition, it should be sufficient to con- 
sider only the static fluctuations of the phase, and we 
may suppress the time-dependence in 9 and retain only 
the spatial fluctuation. The interaction term is then re- 
duced to 

Si = -T J2 mv F (k)-v s (q)^t( k + q, w )Vi(k,a;). 

(2.12) 

Here we take k/m Ri vp(k) where vp(k) is the Fermi 
velocity at the Fermi surface point closest to the wave 
vector k (see Fig. ||), and v s (q) = J d 2 x e™ jq ' x V6>(x)/2m 
is the superfluid velocity. One should bear in mind that 
v s (q) includes both longitudinal and transverse compo- 
nents. 



When the interaction is expressed in form of Eq. (2.12) 
we see that 5*1 is the Doppler shift of the fermion spec- 
trum due to a non-zero superfluid velocity. This becomes 
more explicit if we consider the q — > limit of Si, ne- 
glecting the q dependence of the field operators, so that 
the superfluid velocity can be subsumed into the fermion 
Green's function as 
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where Pf = mvp(k). A semiclassical approximation pf 
this form has recently been used by Franz and MillistZl 
to study phase fluctuations in the underdoped cuprates. 
They incorporated the phase fluctuations by averaging 
the Green's function over a Gaussian distribution of ve- 
locity fluctuations. Our approach is rather different — we 
treat the coupling to the fluctuating phase using a self- 
consistent perturbation theory (see Sec. Ill) in order to 
calculate the fermion self-energy due to phase fluctua- 
tions. 

Finally, we can check if the effective theory correctly 
gives the effective action of the phase fluctuations!^ by in- 
tegrating out the remaining fermionic degrees of freedom. 
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By expanding e Sl to second order in 9, and integrating 
out the fermion fields, we obtain 
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The superfluid density at temperature T i 
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Therefore, for |q| <C A/v F , we recognize that Eq. (2.14) 
has the following form: 
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We see that the stiffness of the phase fluctuations can be 
identified as the superfluid density at T. This result holds 
even when we consider unconventional order parameter 
symmetry. 

Equation (2.16) can serve as a convenient effec- 
tive theory to describe the critical properties .of the 
Berezinskii-Kosterlitz-Thouless (BKT) transition.!^ The 
macroscopic measurable quantities are often insensi- 
tive to the detailed microscopic fermion theory except 
through the parameter n s (T)/m, and the effective de- 
scription of phase fluctuations is enough to characterize 
the transitions in macroscopic quantities. A brief discus- 
sion of the BKT transition and derivation of the effective 
classical action are included in Appendix A for complete- 
ness. We should stress that we are not proposing that the 
superconducting transition in the cuprates is of the BKT 
type; weak interplanar coupling will produce a transition 
in the d — 3 XY universality class, sufficiently close the 
transition. However, outside this regime two-dimensional 
vortex fluctuations will dominate the low-energy physics, 
and we expect the physics discussed in the subsequent 
sections to pertain. Indeed, recent time-domain terahertz 
spectroscopy measurementsE3 of the complex conductiv- 
ity a{oj) of underdoped Bi2Sr2CaCu20s+5, which provide 
a direct probe of the superfluid density, show that the 
dynamics is well described by the BKT picture, lending 
credence to the model developed in this paper. 



To calculate the fermionic self-energy, we have to per- 
form a perturbative expansion in 9 from Eq. (^^), but 
close to T c and for a small superfluid density, there is 
no obvious small parameter. We need a way to select 
a set of meaningful diagrams. The saving grace in this 
case is that the volume (perimeter) of the Fermi surface 
is large, so that we can select the diagrams of leading 
order in l/N(0). One way of keeping track of the dia- 
grams is to introduce n fermion species and perform a 1/n 
expansion and at the end take n — ► 1. This amounts to 
summing over ring diagrams. Other diagrams are smaller 
by factors of X/Uf since only the ring diagrams involve 
summations over the Fermi surface perimeter in this low- 
energy effective theory. In this section we study the effect 
of phase fluctuations on the single-particle properties of 
the fermions using this machinery, both at T < Tbkt 
and T > Tbkt- For the rest of the paper, we rein- 
state the angular dependence of the gap function A(0), 
and perform calculations assuming a d x 2_ y 2--wave sym- 
metry, A(0) = A cos 20, neglecting the temperature de- 
pendence in Ao since it changes smoothly through the 
transition in the underdoped cuprates. 
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FIG. 3. The propagator of the velocity field in the leading 
order in X/Uf. The thin (thick) wiggly line is the bare (full) 
propagator. The ring in the figure is the polarization of the 
fermions. 




FIG. 4. The Dyson equation of fermion Green's function. 
The wiggly line is the velocity propagator. The thin (thick) 
solid line is the bare (full) Green's function. 

We first calculate the velocity- velocity correlation func- 
tion from the effective theory we obtained in Sec. [n] by 
summing over the ring diagrams as in Fig. [| We find 
that 
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III. SINGLE-PARTICLE PROPERTIES OF 
FERMIONS WITH PHASE FLUCTUATIONS 

The underdoped cuprate superconductors can be re- 
garded as doped Mott insulators with a small charge 
carrier density (superfluid density), in other words, 
n s (0)/m ~ xEp where x is the doping concentration. 



xTr[G (k + q,^)G (k,^)]. 

Now we calculate the fermion self-energy correction from 
the Dyson equation, 



S(k,c) = ^[mv F (k)] Q [mv F (k)] 







(3.3) 
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where G is the full fermion Green's function, given self- 



consistently by G 



G,, 



E. As shown in Ap- 



pendix A, the velocity-velocity correlation function can 
be resolved into a longitudinal comp onent C/(q) and a 
transverse component C t (q) [see Eq. ( A14 )]; as a result, 
the self energy can be written as a sum of a contribution 
from longitudinal fluctuations and one from transverse 
fluctuations, 



E(k,w) = E,(k,w)+E t (k,w) 



(3.4) 



where 



£,(k,w) = £[mv F (k) •<7] 2 C ; (q)G(k-q,c), (3.5) 



and 



S t (k, W )=^[mv F (k)xg] 2 C ( (q)G(k-q, W ). (3.6) 



The self-energy has both a momentum and frequency de- 
pendence. Here we focus on the behavior of the self- 
energy near the Fermi surface, assuming that it varies 
smoothly near the Fermi surface, and we neglect the £k~ 
dependence. Therefore the only momentum dependence 
is through the angle cf> on the Fermi surface. 



A. Longitudinal phase fluctuations 

First we study the effect of the longitudinal phase fluc- 
tuations, neglec ting the effect of vortices for now. Using 
Eqs. ( |3.5| ) and (A15), we obtain for the self-energy (see 
Appendix B for details) 
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where uj and A can be calculated self-consistently by 

s(<M- 

(3.8) 



iQ -A((f>) 
-A(</>) iuj 



iuj -A(0) 
-Af» 



By analytically continuing the frequency icj — > uj + irj in 
Eq. (3.S), we can calculate the single-particle density of 
states 
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In order to apply this result to the underdoped 
cuprates, we need to know the temperature dependence 
of the superfluid density n s (T) and the gap A (T). We 
neglect the temperature-dependence of A (T) , since the 



gap smoothly evolves into the pseudogap above T c . A 
reasonable form for the superfluid density that approx- 
imates the BCS theory isJ3 |S (T)/4m rj n s (0)/4m - aT 
where a is order of unity.EIJ At low temperatures, since 
the reduction of the superfluid density is approximately 
linear in temperature, this form is reasonable espe- 
cially for underdoped cuprates where T c is relatively low. 
Throughout this paper we take a = 3/4, but the exact 
value of a is not crucial. Generally the value of a is 
dependent on the value of Ao. 

The curves in Fig. [| show that the peak in the density 
of states (DOS) near ui = Ao is mildly smeared due to the 
phase fluctuations and the effect is more pronounced at 
higher temperatures when the superfluid density is lower. 
Notice that the DOS near uj = is hardly affected even as 
the temperature is raised, and therefore we expect that 
the electronic entropy is only slightly increased due to the 
longitudinal phase fluctuation as illustrated in Sec. IV. 
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FIG. 5. Density of states in the presence of the longitudi- 
nal phase fluctuations at four different temperatures. Here 
t = 4mT/n a (0). Above t — tsKT ~ 0.67, the vortices are 
expected to be more important. Here we take A = Ao/vp. 



B. Transverse phase fluctuations 

Since the two-dimensional superconducting transition 
is driven by vortex-antivortex pair unbinding, we need to 
take into account the effect of vortices, especially above 
the vortex pair unbinding temperature Tbkt^ We as- 
sume that the longitudinal phase fluctuation is less im- 
portant in this temperature regime and consider only 
the effect of the vortex distribution or the transverse 
superfluid velocity (see A ppend ix A). Using the Dyson 



equation similar to Sec. Ill A and neglecting the in- 
dependence in the self-energy, from Eqs. (3.6) and (A16) 
we obtain the following form of the fermion self-energy 
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where q = gj_ + q 2 . Again using Eq. (3.8) we can self- 

consistently calculate ui and A, and subsequently, the 
Green's function when T > Tbkt by analytically contin- 
uing the frequency iui — > lo + ir). 

The ARPES studies of the cuprates show that the 
quasiparticle spectral function broadens dramatically 
whenj-passing from the superconducting to the normal 
state.Q How does vortex unbinding contribute to this 
broadening? In Fig. ^ we show the scattering rate due 
to the vortices, 1/r = Imu) at ui = Aq in the three di- 
rections ((f) = 0, 7r/8, 7r/12). As the temperature is raised 
from Tbkt, 1/t rapidly increases to a value of order 
Ao, which is a rather large quantity. This shows that 
the phase fluctuation is an important ingredient in the 
broadening of the spectral function. As the temperature 
is further raised, we expect that the fluctuations in the 
gap magnitude A will further enhance the spectral width, 
which is beyond the scope of discussion in this paper. 




FIG. 6. The scattering rate (Imu) as a function of 
t = 4mT/n a (0) at u — Ao and in the three directions 
4>\ = 0, <f>2 = 7r/8, 03 = 7r/12, in the presence of unbound 
vortex-antivortex pairs. 

So far we have discussed the spectral properties of the 
neutral fermions. Of much interest is the spectral broad- 
ening of the physical electrons as revealed by the ARPES 
data. For instance, in the absence of the interaction term 
Si in the action, the Green's function for the physical 
electrons factors as 



We may gain insight into the 
spectral function from Eq. (|3 . 1 lh . 



form of the electron 
Above Tbkt, the 

correlation function (e lt ^ x ^ 2 e~ l0{ - ^ 2 ) behaves approx- 
imately-as e~' x '/ 4 ^+, where £+(T) is the BKT correlation 
length,E3 and the electron Green's function is roughly 



<c(x)ct(0)> 



0(x). 



(3.12) 



In the low-energy effective theory we consider here, we 
may think of the quasiparticles as following a quasiclas- 
sical trajectory with a Fermi momentum k^. There- 
fore, the leading contribution to the two-point correla- 
tion function (c(x)c T (0)) comes from the particles with 
the Fermi momentum kj? in parallel with the vector 
±x. Then we may treat x and kp as effectively one- 
dime nsion al£3 The one-dimensional Fourier transform of 
Eq. (3.12) shows that the spectral function of the elec- 
trons is additionally broadened by the width uf/4£ + — 
the electronic spectral width is roughly the sum of 
vf/4£ + and the spectral width of the neutral fermions. 
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FIG. 7. Inverse of BKT correlation length vs reduced tem- 
perature. 

As a result of the rapid spectral broadening above 
Tbkt, the single-particle DOS is significantly modified 
and the singularity at lo = A in the DOS is weakened 
and broadened. Unlike the case of spectral functions, 
the neutral fermion DOS coincides with that of physical 
electrons. Fig. || shows the DOS in the maximum gap 
direction at four different temperatures. When t > 0.67, 
the DOS is rapidly smeared out and the gap is filling in, 
exhibiting a pseudogap-like behavior. 
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FIG. 8. Density of states at <f> = (in the maximum gap 
direction) at four temperatures above Ibkt = 0.67 when the 
unbound vortex-antivortex pairs are taken into account. 



IV. THERMODYNAMIC PROPERTIES 

Experiments on Yo.sCao.2Ba2Cu307_5 show that 
there is a BCSJikc discontinuity in the specific heat at op- 
timal doping.tl However, in the underdoped regime, the 
superconducting transitions are marked by weaker and 
broader peaks in the specific heat. With this puzzling 
feature in mind, we consider the effect of the phase fluc- 
tuations on specific heat in this section. 

Before we proceed, we consider one consequence of 
the pseudogap phenomenon in underdoped cuprates. We 
can estimate the specific heat change in the presepce of 
the pseudo gap using the BCS mean-field theory;EJ the 
superconducting-normal free energy difference is 
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Here g is the BCS pair coupling and A s and A are the 
superconducting gag and the total electronic gap respec- 
tively. The ARPESQ and STSE3 data show that the elec- 
tronic gap evolves smoothly through T c , from the super- 
conducting gap to the pseudo gap. Since A(T) does not 
rapidly change at T c and Eq. (4.2) is expressed in terms 
of A(T) only, d(l/\g\)/dA s is a very small number and 
there is no dramatic increase in entropy at T c . Therefore, 
we do not expect a BCS-like specific heat jump. 

From the above discussion we see that most of the 
characteristic changes in the specific heat at T c should 
come from the order parameter fluctuations. Here we 



consider the effect of the phase fluctuations at both T > 
Tbkt and T < Tbkt- First we consider the electronic 
entropy change due to the longitudinal phase fluctuations 
which is more important at T < Tbkt- From the Green's 



function and the density of states we obtained in Eq. (3.7) 
and Eq. (3J3), we may approximately express the entropy 
as 



+ [1-/M]ln[l-/M]}, (4.3) 

where N(uj) is the single-particle density of states and 
f(x) = l/(e^ x + 1). We can compare the entropy thus 
obtained with the mean-field BCS entropy, as illustrated 
in Fig. [| We observe that the electronic entropy change 
is so minute that it may be hardly visible in the specific 
heat. The small entropy increase is due to the fact that 
the density of states is affected by the phase fluctuations 
only when iv « Ao (see Fig. |J); the electronic entropy is 
insensitive to the longitudinal phase fluctuations. 




FIG. 9. The reduced electronic entropy a = ( ^"ff ) & vs 
the temperature, a has the meaning of the entropy in an 
area of 1/kp normalized by the maximum gap energy Ao. 
The dotted line is the electronic entropy obtained using the 
d-wave BCS theory, and the solid line is the entropy in the 
presence of the longitudinal phase fluctuations. Here the case 
of n 3 (0)/4m = Ao/3 is illustrated. 

Now we turn to the behavior of specific heat in the 
presence of unbound vortex-antivortex pairs at T > 
Tbkt- In the 2D XY model, the vortex-antivortex dis- 
sociations cause an unobservable singularity at the tran- 
sition temperature, followecL-by a broad bump at a tem- 
perature higher than Tbkt^ In conventional supercon- 
ductors, the bump in the specific heat due to the vortex 
unbinding may be unobservable — the density of disso- 
ciated vortices can be only so large before they satu- 
rate the whole sample since the vortex core size is much 
larger than the electronic mean spacing. However, in the 
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cuprate superconductors, the vortex core size is only a 
few times the lattice spacing, so the vortex dissociation 
may give a sizable contribution to the specific heat. To 
reduce the theory down to the effective vortex degrees of 
freedom, we integrate out the fermion fields and express 
the action in the form in Eq. ( A7 ) . Then we can evaluate 
the specific heat by 



C 



kB ±( T *± T 
dT\ dT 

where the reduced free energy is defined as 



T = In Z = In 



Dn e 



(4.4) 



(4.5) 



The free energy is a function of the temperature through 
K(0) and y(0) as defined in Appendix A. We evahwte f 
following the prescription by Berker and Nelson,c3 who 
computed the specific heat in the 2D XY model. The 
only difference here is that in our model the superfluid 
density and the vortex core energy have a temperature 
dependence whereas in the spin XY model the spin stiff- 
ness JT is a constant. 




FIG. 10. The reduced heat capacity per area of a vortex 
core (c = (^c/^s) C) vs temperature (t = 4mT/n s (0)). There 
is a very weak singularity at t = 0.67. 

The result is shown in Fig. [l^. The heat capacity per 
vortex core area ^ is peaked at t « 0.8 and the peak 
height about 13fcs. Now we may cast this result on the 
electronic length scale. Since fcpfic ~ 10 in cuprates, the 
heat capacity per unit cell of area 1/kp is about O.lfcs, 
which is a large number compared to the electronic spe- 
cific heat. In conventional low-T c superconductors, there 
is a factor of 10~ 6 reduction in the peak since kpt;^ 10 3 . 
The specific heat peak may be practically unobservable 
in this case. We conclude that if there is a vortex pair un- 
binding transition in underdoped cuprates, its effect on 
thermodynamic quantities would be visible, although the 



actual specific heat measurement may be complicated by 
other sources of thermal fluctuations. 



V. CONCLUSIONS 

The low-energy effective theory which we have applied 
to the underdoped cuprates is a convenient model which 
describes a complex many-body system (whose micro- 
scopic details are unknown) in terms of a few macroscopic 
parameters such as the superfluid density, effective mass, 
and the magnitude of the superconducting gap. Utiliz- 
ing this phenomenological theory, we are able to calculate 
the spectral properties in the presence of phase fluctua- 
tions without regard to the microscopic details. By as- 
suming a robust quasiparticle gap and a small superfluid 
density, as appropriate for the underdoped cuprates, we 
found that the vortex pair unbinding transition causes 
pronounced effects on the single-particle properties: (i) 
The peak in the spectral function is dramatically broad- 
ened as the temperature is raised through Tbkt] (ii) The 
density of states is widely smeared out at T > Tbkt, 
significantly deviating from the usual superconducting 
state although the gap magnitude is robust. These two 
features are reminiscent of the pseudogap behavior in uru 
derdoped cuprates, as observed in the ARPES0 and STScI 
experiments. It is remarkable that we can reproduce a 
few of the pseudogap-like properties by including only 
the phase degree of freedom. This indicates that along 
with the fluctuations in the order parameter amplitude, 
the phase fluctuations may be an essential ingredient in 
the pseudogap phenomena. For instance, phase fluctua- 
tions may better explain the fact that the superconduct- 
ing gap nodes evolve into extended gapless Fermi arcs 
in the pseudogap phase since the phase fluctuations are 
more significant near the nodes (Fig. ^). The amplitude, 
fluctuations alone may not be sufficient to explain this.E£l 

The role of phase fluctuations in the pseudogap regime 
has recently been explored in several approaches which 
bear some similarity to ours. As mentioned in Sec. II, 
Franz and MillisEZI have also studied the effect of classical 
phase fluctuations on the spectroscopic properties of the 
underdoped cuprates, by coupling the d-wave quasipar- 
ticles to the fluctuating supercurrents due to unbound 
vortex-antivortex pairs. While our approach to treat- 
ing this coupling is quite different from the approach of 
Franz and Millis, our conclusions are similar — that trans- 
verse phase fluctuations are important in determining 
the spectral properties aboiie T c . On a different tack, 
Balcnts, Fisher, and NayakEZI have recently proposed a 
zero-temperature "nodal liquid theory" in which there is 
a transition between the superconducting phase and the 
nodal liquid (tentatively identified as the ground state of 
the pseudogap regime) which is driven by quantum fluc- 
tuations of the order-parameter phase. The properties 
of this model at non-zero temperatures have yet to be 
investigated. 
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Before closing we should note that our interpretation 
of the pseudogap regime as the remnant of a supercon- 
ducting phase which has been destroyed by phase .fljjctu- 
ations is got without criticism. Geshkenbein et alEj and 
RanderiaB have argued that the fluctuation effects should 
be significant over a much wider temperature range than 
is observed experimentally, and that strong pairing cor- 
relations (beyond BCS theory) need to be incorporated 
into any model of the pseudogap regime. A completely 
different explanation of therpseudogap behavior has been 
proposed by Lee and Wcn,Eil who argue that the super- 
conducting phase is destroyed by the thermal excitation 
of nodal quasiparticlcs, which quenches the superfluid 
density while retaining the quasiparticle gap. Finally, 
we have neglected the Coulomb interaction in our model; 
including it would presumably suppress the longitudinal 
phase fluctuations, while leaving the transverse fluctua- 
tions unaffected. Ell All of these criticisms merit further 
investigation in a more detailed study of the pseudogap 
phenomenon. 
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APPENDIX A: LONGITUDINAL AND 
TRANSVERSE FLUCTUATIONS OF THE 
SUPERFLUID VELOCITY 

For the sake of completeness, in this Appendix we 
have collected together some of the important results 
regarding.— .the Berezinskii-Kosterlitz-Thouless (BKT) 
transitionlij which pertain to our calculation of the 
fcrmion self-energy. A more detailed discussion can be 
found in Ref. First, we resolve the phase field into 
longitudinal (spin wave) and transverse (vortex) com- 
ponents. This is done by writing 8 = 6 a + 9 V , with 
9 a the analytic spin-wave contribution and 6 V the sin- 
gular vortex contribution. The superfluid velocity can 
then be decomposed as v s = (l/2m)V6* = v' + v*, 
with = (l/2m)V0 a and v* = (l/2m)V0„, such that 
Vxv^=0 and V 
then becomes 



v* = 0. The classical phase action 



mn s (T) 

2T 
mn s (T) 

2T 



j2 2 

a xv a 



since the cross term is zero as shown below. 

The curl of the transverse component is the total 
vorticity; for singly quantized vortices with circulation 
±27r?i/2m, this can be written as 



V x v s = — n (x)z, 



where the vortex density n(x) is 

ra(x) = y^(5(x-Xj), 



(A2) 



(A3) 



with vortices at positions {x^ } and with vortex "charges" 
qi = ±1 (here we consider only neutral vortex configu- 
rations such that = 0)- Using the fact that v* is 
transverse, 

/9ir 
d 2 x'G(x-x')—n(x')z, (A4) 

where V 2 G ~ — <5(x — x'); Fourier transforming, 

vS(q) = n(q). A5 

The analogous expression for the longitudinal component 
is 



Vi(q) = ^.(q). 



(A6) 



It is e asy to verify that the cross term v*(q) • v' (q) in 
Eq. (Al) is zero. The final expression for the vortex 
contribution to t he c lassical ph ase action is obtained by 
substituting Eq. ( |A5| ) into Eq. flAlj ): 



S v =irK(Q) J J d 2 a;dVn(x)n(x')ln(|x-x'|/£c) 

|x-x'|>£ c 

+ lny(0) J d 2 xn 2 (x) 



(f 



-8^ c 2 ln[y(0)] n(qM-q), (A7) 



where K(0) — n s (T)/AmT is the bare spin stiffness, 
?/(0) = exp(— Ec/A-kT) is the bare fugacity with £ c = 
TK(0)\uk the vortex core energy, £ c ~ vf/Aq is the 
vortex core size, and n is the ratio of the magnetic pen- 
etration depth to the coherence length. 



t\2] 



d 2 x{(v i s y + (vi) 



(Al) 
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y(i) 




FIG. 11. The RG flow of 2/(0)- The dashed 

curve is the locus of initial values (K~ (0),y(0)). The 
asterisk (*) marks the vortex pair unbinding transition 
temperature (Tbkt) which in this model corresponds to 
4mTBKT/n a (0) ~ 0.67. Above this transition temperature, 
(if _1 (0), y(0)) rides on the outgoing curves and (K -1 (l),y(l)) 
values diverge in the I — > oo limit. 



phase action, Eq. (A7), is quadratic in the vortex den- 
sity, the correlation function is 



,-2Z* 



T 2 (q, K(0),y(0)) 



4ir 2 K(l*)/(e l *q) 2 - 8tt£ 2 \ny{l 

If we choose I* such that 

4 7 r 2 J fY(/*)e c 2 = -87re c 2 lny(r), (All) 

then the correlation function takes the particularly sim- 
ple form 



(Al 



T 2 (q,K(0),y(0)) = 



1 



1 



4ir 2 K(l*) l/q 2 + £ 2 e 21 ' 



(A12) 



Both I* and K(l*) arc determined f rom a numerical in- 
tegration of the RG equations, Eq. (AS). 

Using these results we can calculate the velocity- 
velocity correlation function for T > Tbkt, 



«(q)^(-q))=Q(q)g a g /3 

+C t (q) (S a0 - q a qp), 

where the longitudinal component is 



(A13) 



In this model there is a critical temperature Tbkt 
above which the two-dimensional quasi-long-range order 
is destroyed by thermal fluctuations. At T < Tbkt, 
the vortices are bound into neutral pairs (when there is 
no external magnetic field), and the superfluid density 
is non-zero. In this temperature regime, consideration 
of the soft longitudinal phase fluctuation a is sufficient. 
When T > Tbkt, the bound pairs of vortices dissociate 
and free vortices proliferate as the temperature is raised, 
driving the superfluid density to zero. This picture is 
made quantitative using the RG flow equationst^l for the 
coupling constants (K(l),y(l)) at the length scale e l (see 
Fig. Ill): 



dK-\l) 
dl 

dy(l) 
dl 



= 47r 3 y 2 (0, 

= [2-irK(l)]y(l). 



(A8) 



To calculate the fermion self-energy we need the 
vortex-vortex correlation function (n(q)n(— q)) = 
r 2 (q, K(0), ?/(0)), which can be calculated approximately 
using the Debye-Huckel theory developed by Halperin 
and Nelson. t£l To implement this, we first note that un- 
der scale transformations the correlation function scales 



r 2 ( q ,ir(o),y(o)) 



,-2i 



T 2 (e l q,K(l),y(l)). (A9) 



Now we integrate the RG flow equations out to a length 
scale I* where the density of free vortices is high; at 
this scale we may integrate over the coarse-grained vor- 
tex density rather than the individual vortex coordinates 
(the Debye-Hiickel approximation). Since the classical 



Ci(q) = ^<0«(q)0a(-q)> 
T 



(A14) 



Am 2 n s (T) ' 
and the transverse component is 

7T 2 (n(q)n(-q)) 



Q(q) 



q* 

1 



(A15) 



1 



Am 2 K{l*) l + £, 2 e 2l 'q 2 ' 



APPENDIX B: DETAILS OF THE 
CALCULATION OF THE SELF-ENERGY 

In this Appendix we provide some of the details of the 



calculation of the self energy in Eqs. (p.7[) and (3 



Near the Fermi surface, neglecting the ^-dependence 
in E(k, w), 



iui - £ k -A(0) 
-A((p) iw + ^k 



(Bl) 



where <fi is the angle on the Fermi surface. From Eq. 
(p|), when |k| = k F , 



iCj + <S£ q A(<p) 



E,(^«) « -J2 [mv F (k) -qfciiq) 



+ (<5£ q ) 2 + A 2 (cb) ' 



A(</>) iu - <f£ q 
(B 
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where <5£ q = Vf(k) • q. Now, using the integral 

c 2 1 



dx dy 



x 2 + y 2 x 2 + a 2 



X /* oo 2 

X 



1 



x J-oo x +V x + a 



7T In 



(B3) 



we obtain Eq. (3.7) 
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